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Abstract 

(N ■ 

!> ■ The maximum overproduction of helium-4 in cosmological nucle- 

osynthesis with active-sterile neutrino oscillations, v e <-> u s , efficient 
\ after decoupling of electron neutrino, is analyzed. The kinetic effects 

ON ' on primordial nucleosynthesis due to neutrino spectrum distortion, 

caused by oscillations, are precisely taken into account. 
\ The maximum overproduction of primordial 'Tie as a function of oscil- 

lation parameters is obtained from the analysis of the kinetics of the 
nucleons and the oscillating neutrinos, for the full range of parame- 
q , ters of the discussed oscillation model. A maximum relative increase 

\ of ^e, up to 14% for non-resonant oscillations and up to 32% for 

resonant ones is registered. Cosmological constraints on oscillation 
parameters are also discussed. 
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1 Introduction 



In this work I study the maximum overproduction of Tie in Big Bang Nucle- 
osynthesis (BBN) with electron-sterile neutrino oscillations v e «-> u s . 

The positive indications for oscillations, obtained by the neutrino exper- 
iments (SuperKamiokande, SNO, Soudan 2, LSND, etc.) turned the subject 
of neutrino oscillations into one of the hottest points of astrophysics and neu- 
trino physics. The solar neutrino problem, the atmospheric neutrino anomaly 
and the positive results of LSND experiments may be naturally resolved by 
the phenomenon of neutrino oscillations, implying nonzero neutrino mass 
and mixing. Although sterile neutrino impact in oscillations explaining at- 
mospheric and solar neutrino anomalies, is strongly constrained from the 
analysis of experimental oscillations data, still some small fraction of v s may 
participate in oscillations. Hence, it is interesting to study cosmological ef- 
fects of such oscillations. 

Cosmological nucleosynthesis with neutrino oscillations was studied in 
numerous publications, discussed in detail in ^ . In these publications the 
central goal was to obtain cosmological constraints on oscillation parameters. 

The overproduction of Tie itself was not studied in detail until now. Such 
a study may be of interest for constructing of alternative cosmological mod- 
els, for constraining galactic chemical evolution, it may be useful for non- 
standard models predicting active-sterile oscillations, like models with extra 
dimensions, mirror world particles, etc. 

Here I address that question of the possible maximal overproduction 
of Tie due to oscillations. For that purpose the case of v e <-> u s , oscilla- 
tions, effective after the electron neutrino decoupling from the plasma (i.e 
for Sm 2 sin 2 2$ < 10 _7 eV 2 ), is the most suitable. In that case due to the fact 
that sterile neutrino state is usually less populated than the active one at 
the start of oscillations, the oscillations may cause strong distortion of the 
electron neutrino spectrum, which affects the kinetics of the nucleons freezing 
before nucleosynthesis, and hence, the primordial production of Tie. And, 
as will be shown by the numerical analysis, this kinetic effect of oscillations, 
may be much bigger than the one corresponding to an increase in the total 
energy density due to an additional neutrino flavor 5N p , mainly considered 
in literature. 

We analyze Tie primordial production, taking into account all known ef- 
fects of v e <-> v s oscillations on the primordial synthesis of Tie. The produc- 
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tion of ^e was calculated in the non-resonant and resonant oscillation cases. 
In both cases strong overproduction of ^e was found possible - up to 14% 
and 32%, correspondingly. Thus in the discussed model of non-equilibrium 
oscillations the maximum overproduction of ^e corresponds to an increase 
of the neutrino effective degrees of freedom SN^ X ~ 6. 

The oscillation effects on BBN and their description are briefly reviewed 
in the next section. The kinetic approach, the results on 4 He primordially 
produced abundance in the presence of oscillations, and the cosmological 
constraints on oscillation parameters are discussed in the last section. 

2 Cosmological Nucleosynthesis with Neutrino 
Oscillations 

2.1 Standard Cosmological Nucleosynthesis 

According to the standard BBN, during the early hot and dense epoch of the 
Universe, the light elements D, ^e, ^e, \ii were synthesized successfully. 
The most reliable and abundant data are now available for Tie. This fact 
and its relatively simple chemical evolution make ^e the preferred element 
for the analysis of the oscillations effect on BBN. 

The contemporary values for the mass fraction of Tie, Y p , inferred from 
observational data, are in the range 0.238-0.245 (the systematic errors are 
supposed to be around 0.007) |§. 

Primordially produced ^e abundance Y p , is calculated with great preci- 
sion within the standard BBN model [|J. According to the standard BBN 
'Hie is a result of a complex network of nuclear reactions, proceeding after 
the neutron-to-protons freezing. It essentially depends on the freezing ratio 
(n/p)f. The latter is a result of the freezing of the weak processes: 

v e + n <-> p + e~ 

e + + n <->• p + z> e , (1) 

which maintained the equilibrium of nucleons at high temperature (T > 1 
MeV). Their freeze-out occurred when in the process of Universe cooling 
the rates of these weak processes, T w , became comparable and less than the 
expansion rate H(t): 
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r w ~ G 2 F E 2 V N V < H(t) ~ v^77 t2 

Thus, the produced Tie is a strong function of the number of the effective 
degrees of freedom at BBN epoch, g e jf = 10.75 + 7/45N p . 

Y p depends also on the electron neutrino spectrum and number density, 
and on the neutrino-antineutrino asymmetry, which enter through T w . In 
the standard BBN model three neutrino flavors, zero lepton asymmetry and 
equilibrium neutrino number densities and spectrum distribution are postu- 
lated: 

n^(E) = (l + eME/T))- 1 . 

Almost all neutrons, present at the beginning of nuclear reactions, are 
sucked into Tie. So, the primordially produced mass fraction of Tie can be 
approximated by 

Y p ~ 2(n/p) f /(l + n/p)/exp(-t/r n ). 

where r n is the neutron mean lifetime. 

The theoretical uncertainty of the precisely calculated Y p is less than 0.1% 
{\SY p \ < 0.0002) within a wide range of rj. So, in the standard BBN scenario, 
where r n , and g e ff are fixed, Y p , as well as the rest cosmo logically produced 
light elements are functions of only one parameter - the baryon-to-photon 
ratio rj. 

Deuterium measurements in pristine environments towards low metal- 
licity quasar absorption systems at very high z ~ 3 provide us with the 
most precision determination of the baryon density ||, giving the value: 
rj = (5.6 ± 0.5) x 10~ 10 . Recently, the baryon density was also determined 
from observations of the anisotropy of the cosmic microwave background 
(CMB) by DASI [§, BOOMERANG §, g, CMB § and MAXIMA exper- 
iments PP|| . For the combined analysis of the data see also ]TTj, [T2]]. The 
CMB anisotropy data is in remarkable agreement with the baryon density 
determined from deuterium measurements and BBN. 

So, the predicted primordial Tie abundance Y p is in accordance with the 
observational data and is consistent with deuterium measurements. 
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2.2 Effects of neutrino oscillations on nucleosynthesis 

In case neutrino oscillations are present in the Universe primordial plasma, 
they lead to changes in the Big Bang Nucleosynthesis, depending on the 
oscillation channels and the way they proceed. The effect of flavor neutrino 
oscillations on BBN is negligibly small because the temperatures and hence, 
the densities of the neutrinos with different flavors are almost equal. 

Active-sterile oscillations, however, are capable to shift neutrino num- 
ber densities and spectrum from their equilibrium values. Besides, they 
may change neutrino-antineutrino asymmetry and excite additional neutrino 
types. Thus, the presence of neutrino oscillations invalidates the main BBN 
assumptions about three neutrino flavors, zero lepton asymmetry and equi- 
librium neutrino energy distribution. 

Qualitatively, the oscillations effects on the nucleosynthesis of 'Tie may 
be illustrated as follows: 



• excitation of additional degrees of freedom 



It is known that active-sterile oscillations may keep sterile neutrinos in ther- 
mal equilibrium fl3| or bring them into equilibrium in case they have already 



decoupled. The presence of light steriles in equilibrium leads to an increase 
of the effective degrees of freedom during BBN and to faster Universe ex- 
pansion H(t) ~ g l Jj 2 f and earlier n/p- freezing, Tf ~ (gef/) 1 / 6 , at times when 
neutrons were more abundant |T3|, \lT 



n/p~ exp [-(m n - m p )/T f ] 

This effect leads to 5% "tie overproduction corresponding to one additional 
neutrino type brought into equilibrium by oscillations. 



• distortion of the neutrino spectrum 



Much stronger effect of oscillations may be achieved in case of oscillations 
between initially empty sterile neutrino state and electron neutrino. The 
non-equilibrium initial condition leads to spectrum distortion of the active 
neutrino due to oscillations. 

Since oscillation rate is energy dependent T ~ 5m 2 /E the low energy 
neutrinos start to oscillate first, and later the oscillations become noticeable 
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for the more energetic neutrinos. Due to that, the neutrino energy spectrum 
n v {E) may strongly deviate from its equilibrium form [15|. This spectrum 
distortion affects the kinetics of nucleons freezing - it leads to an earlier 
n/p-freezing and an overproduction of Tie yield. 

The effect can be easily understood having in mind that the distortion 
leads both to a depletion of the active neutrino number densities N u 



16 



&EE 2 n v (E) 



and to a decrease of the mean neutrino energy, which reflects into a decrease 
of the weak rates r w ~ E 2 N U , and hence, to an overproduction of Tie pri- 
mordial abundance. 

The decrease of the electron neutrino energy due to oscillations into 
low temperature sterile neutrinos, has also an additional effect: Due to the 
threshold of the reaction converting protons into neutrons, when neutrinos 
have lower energy, protons are preferably produced in reactions (1), which 
may lead to an underproduction of Tie [[I?]. However, this turns to be a 
minor effect. 

The effect of spectrum distortion was numerically analyzed for hundreds 
combinations of mass differences and mixing angles in previous studies of 
active-sterile oscillations [15, [Tj3, 19, 2C]. It was proved important both in 
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and in the non-resonant oscillations case 
In the discussed here scenario it is the dominant effect and leads to over- 
production of the primordial Tie abundance. 



neutrino-antineutrino asymmetry growth 



The idea of neutrino-antineutrino asymmetry generation during the resonant 
transfer of neutrinos was first proposed in ref. [f21|. Dynamically produced 



asymmetry exerts back effect to oscillating neutrino and may change its oscil- 
lation pattern p2l E3j , it may suppresses oscillations at small mixing angles, 



leading to weakening the oscillations effect on BBN, i.e. less overproduction 
of -H< • 

The effect of the oscillations generated asymmetry on Tie for the discussed 
here model was analyzed for hundreds of 5m 2 — •& combinations in 
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2 There were independent studies of the asymmetry growth in active-sterile oscillations 
for the case of large mass differences along the lines of the pioneer work [E4| . 
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3 Helium-4 overproduction due to v e <-> v s 
neutrino oscillations 



We have provided an exact study of all the oscillation effects on the primordial 
production of 4 He. 

3.1 The required kinetic approach 

For the analysis of the non-equilibrium picture of active-sterile neutrino os- 
cillations, producing non-equilibrium neutrino number densities, distorting 
neutrino spectrum and generating neutrino-antineutrino asymmetry a self- 
consistent numerical analysis of the evolution of the nucleons and the oscil- 
lating neutrinos in the high temperature Universe was provided. 

Exact kinetic equations for the nucleons and for the neutrino density ma- 



trix in momentum space |18| were used. This allowed to describe precisely the 
kinetic effects of oscillations on helium production due to spectrum distortion 
at each neutrino momentum. 

The equation for the neutron number densities in momentum space n n 
reads: 

(dn n /dt) = Hp n (dn n /dp n ) + 

+ J dtl(e~,p, u)\A(e~p -> vn)\ 2 [n e -n p {l - Pll)~ 

-n n pLL(l - n e -)\ 

-J dtt(e + ,p, v)\A{e + n -> pv)\ 2 [n e +n n (l - Pll)~ 
-n p p LL (l - n e +)] . (2) 

where dfl(i, j, k) is a phase space factor and A is the amplitude of the cor- 
responding process, neutrino pn and antineutrino number densities pn at 
each integration step of eq. (2) are taken from the simultaneously performed 
integration of the set of equations for neutrino density matrix (see ref. ||18|| ). 

The equation provides a simultaneous account of the different compet- 
ing processes, namely: neutrino oscillations (entering through pll and Pll), 
Hubble expansion (first term) and weak interaction processes (next terms). 

The numerical analysis was performed for the temperature interval [2 MeV, 
0.3 MeV]. 
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We have found that for a wide range of oscillation parameters the spec- 
trum distortion of electron neutrino is considerable during the period of nu- 
cleons freezing. Hence, usually the kinetic effect of oscillations due to electron 
neutrino energy spectrum distortion plays the dominant role in the overpro- 
duction of helium. In Fig. 1 the evolution of the energy spectrum of the 
electron neutrino through the period of nucleons freezing is illustrated. In 
the Figs, la, lb and lc the energy spectra at different characteristic tempera- 
tures, correspondingly T = 1 MeV, T = 0.7 MeV and T = 0.5 are presented. 
By the dashed curve the equilibrium spectrum at the given temperature is 
given for comparison. The spectra are calculated for oscillation parameters 
5m 2 = 10~ 7 eV 2 and sin 2 2$ = 0.1. 

It is seen that oscillations affect non trivially the neutrino spectrum and 
distort strongly its equilibrium form. So, it is not possible to describe cor- 
rectly the spectrum distortion due to oscillations simply by shifting the effec- 
tive temperature of the neutrino and accounting only for the depletion of its 
total number density, considering its spectrum equilibrium as in refs. [^, pTf] . 

For the proper description of the spectrum distortion in the case of non- 
equilibrium electron-sterile oscillations, studied here, the evolution of the 
neutrino ensembles should be explored using neutrino density matrix in mo- 
mentum space, allowing to describe the evolution of neutrino at each mo- 
mentum. 

The neutrino-antineutrino asymmetry in the resonant case grows up to 
5 orders of magnitude from its initial value (taken to be of the order of the 
baryon asymmetry). So, this asymmetry influences BBN only indirectly - 
through oscillations. This oscillations generated asymmetry leads to a de- 
crease in the produced Tie compared with the case without asymmetry ac- 
count. However, its effect comprises only up to about a 10% of the total 
effect. 



3.2 Maximum helium-4 overproduction 

The overproduction of the primordial Tie, SY P = Y° sc — Y p in the presence of 
v e *-* v s oscillations was calculated for mass differences 5m? < 10~ 7 eV 2 and 
< $ < 7r/4. Several hundreds of Sm 2 — d combinations were explored. 

We have used the data from the precise calculations of the n/p-freezing 
provided for the non- resonant case in pOf and for the resonant case in [|i~9|1 . 



As far as it is the essential for the production of Tie. The neutron decay was 
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accounted adiabatically till the beginning of nuclear reactions at about 0.09 
MeV. 

We have found that the effect of oscillations becomes very small (less 
than 1%) for small mixings: as small as sin 2 2$ = 0.1 for 5m 2 = 10~ 7 eV 2 , 
and for small mass differences: 5m 2 < 10~ 10 eV 2 at maximal mixing. For 
very small mass differences 5m 2 < 10~ u eV 2 , or at very small mixing angles 
sin 2 2t? < 10 -3 , the effect on nucleosynthesis becomes negligible. 

In the non-resonant case the oscillation effect increases with the increase 
of the oscillation parameters, hence it is maximal at maximal mixing and 
greatest mass differences. In Fig. 2 (the lower curve) the maximal relative 
increase in the primordial "Tie as a function of neutrino mass differences at 
maximal mixing: 5Y™ ax /Y p = 5Y° SC /Y p (5m 2 )\g =7T /4 : is presented. 

It is seen that for maximal mixing, the oscillation effect becomes greater 
than 5% (the one corresponding to one additional neutrino type) already at 
\5m 2 \ > 3 x 10~ 9 eV 2 (in the resonant case) and 5m 2 > 6 x 10~ 9 eV 2 (in 
the non- resonant one). It continues to grow up with the increase of the mass 
differences, and at \5m 2 \ ~ 10 -7 eV 2 is several times bigger: 5Nj™% x ~ 3 in 
the non-resonant case and SN^ X ~ 6 for the resonant one. 

Further increase of the mass differences, however, will lead to oscillations 
effective before electron neutrino decoupling, and therefore, to a smaller spec- 
trum distortion effect, because the interactions with the plasma will lead to 
faster thermalization of the sterile state. Hence, the effect on helium will 
decrease with further increase of |<5m 2 | and finally reach an overproduction 
of 5% again, corresponding to a full thermalization of the sterile state and 
its equilibrium spectrum. 

In Fig. 3 we present a combined plot (for the resonant and the non- 
resonant oscillation case) of 5Y P dependence on the mixing angle for 5m 2 = 
10~ 7 eV 2 and 5m 2 = 10~ 8 eV 2 . While in the non-resonant case the oscillations 
effect increases with the increase of the mixing (see l.h.s of Fig. 3.), in the 
resonant case for a given Sm 2 there exists some resonant mixing angle, at 
which the oscillations are enhanced by the medium (due to the MSW effect), 
and hence, the overproduction of ^He is greater than that corresponding to 
the vacuum maximal mixing angle. This behavior of the helium production 
on the mixing angle is illustrated in the r.h.s. of the figure. 

The upper curve in Fig. 2 shows the maximal relative increase in the 
resonant oscillations function of mass differences. Each maximum 

^He value corresponds to the resonant mixing angle for the concrete mass 
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difference: 5Y p max /Y p = Y° sc (6m 2 , As can be seen from Figs. 2 and 3, 

a considerable overproduction can be achieved: in the resonant case up to 
32% and in the non-resonant one - up to 14%. So, the net effect of spectrum 
distortion of electron neutrino due to oscillations on the production of Tie 
may be considerable and several times larger than the effect due to excitation 
of one additional neutrino type. 

Several words are due to cosmological constraints on neutrino oscillation 
parameters, calculated in the discussed model of oscillations. Due to the 
strong kinetic effect of neutrino spectrum distortion caused by active-sterile 
oscillations discussed, the obtained constraints in that model of oscillations 
are more stringent than those obtained in pioneer works, not accounting for 
the kinetic effects of oscillations. Hence, according to these more precise 
studies of spectrum distortion effects of oscillations, the 5Y p /Y p < 3% limit 
excludes almost completely the LOW electron-sterile solution to the solar 
neutrino problem [^5j, in addition to the excluded sterile LMA solution in 
previous investigations. I.e. it is more than an order of magnitude more 
restrictive to the mass differences. 

This study has shown that considerable Y p overproduction may result 
from the electron neutrino spectrum distortion due to u e «-> v s oscillations. 
The overproduction is maximal for the case of initially empty sterile neu- 
trino state, considered here. The dependence of Tie overproduction on the 
degree of population of the sterile neutrino state before v e <-> v s oscillations 



is considered in a separate paper pgj . Bigger initial population of v s leads to 
a smaller spectrum distortion of u e , and hence smaller kinetic effect on pri- 
mordial nucleosynthesis 5Nki n and higher energy density due to the increase 
of the number of degrees of freedom 5N P . The interplay between the two 
effects, however, is such that the overproduction of Tie is smaller than in the 
case when initially 5N p = 0, discussed here. 



4 Conclusions 

The primordial production of Tie in the presence of v e <-» v s oscillations, effec- 
tive after electron neutrino decoupling was analyzed. A precise quantitative 
study of the maximum overproduction of Tie, accounting for all oscillations 
effects is provided. It was shown that the considerable spectrum distortion 
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of the electron neutrino caused by oscillations, which effects the kinetics of 
the neutron-proton transitions during nucleons freezing, plays the dominant 
role in the overproduction of 4 tie. 

Enormous overproduction of ^e (up to 32% in the resonant case and 
14% in the non-resonant case) was found possible in case the sterile neutrino 
state was empty at the start of oscillations. 

The results of this analysis can be useful for constraining nonstandard 
physics, predicting active-sterile neutrino oscillations, like extra- dimensions, 
producing oscillations, supernova bursts employing oscillations, etc.. It can 
be of interest also for models of galactic chemical evolution. 
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Figure la: The figure illustrates the degree of distortion of the electron 
neutrino energy spectrum x 2 pll(x), where x = E/T, caused by oscillations 
with mass difference |<5m 2 | = 10~ 7 eV 2 and sin 2 2$ = 0.1 at a character- 
istic temperature 1 MeV. The dashed curve gives the equilibrium neutrino 
spectrum for comparison. 
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Figure lb: The figure illustrates the degree of distortion of the electron 
neutrino energy spectrum x 2 pll(x), where x = E/T, caused by oscillations 
with mass difference \5m 2 \ = 10~ 7 eV 2 and sin 2 2i? = 0.1 at a characteristic 
temperature 0.7 MeV. 
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Figure lc: The figure illustrates the degree of spectrum distortion of 
the electron neutrino caused by oscillations with mass difference Sm 2 = 10 -7 
eV 2 and sin 2 2$ = 0.1 at a characteristic temperature 0.5 MeV. 
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Figure 2: Maximum primordial "tie abundance for the resonant (upper 
curve) and the non- resonant oscillation case (lower curve), as a function of the 
neutrino mass differences. The non-resonant case is calculated at maximum 
mixing, while in the resonant case the helium abundance is calculated at the 
resonant mixing angle for the corresponding mass difference. 
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Figure 3: The dependence of the relative increase of primordial helium 
on the mixing angle for the resonant (r.h.s.) and non-resonant (l.h.s.) oscil- 
lation case. The upper curve corresponds to 5m 2 = 1CT 7 eV 2 , the lower one 
to 5m 2 = 10~ 8 eV 2 . 
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